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Motivated by the large interest in the non-equilibrium dynamics of low-dimensional quantum 
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which is based on time-dependent current-density-functional theory, is applicable well beyond the 
linear-response regime and produces both spin-charge separation and spin-drag-induced broadening 
of the spin packets. 
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Introduction — The dynamics of low-dimensional quan- 
tum many-body systems is a topic of great theoretical 
and experimental interest. With the recently acquired 
capability to store and control ultracold atomic gases [T], 
it has become possible to study in the laboratory some 
ideal systems which were once accessible only to theo- 
retical investigation. In particular, it is now possible to 
confine atoms along a waveguide, thus realizing for the 
first time the elusive one-dimensional Fermi and Bose 
gases [2]. 

One-dimensional (ID) Fermi systems are theoretically 
interesting because they present us with a breakdown of 
the Landau Fermi liquid paradigm [3J. The low-energy 
excitations of such systems are collective, while those of a 
conventional Fermi liquid are single-particle like; further- 
more there is complete separation between charge and 
spin excitations, which in an ordinary Fermi liquid are 
tied together in the "quasiparticle" . 

Fully-microscopic calculations of collective spin and 
charge dynamics in ID Fermi systems have recently 
been performed by Kollath et al. [A by means of the 
time-dependent density-matrix renormalization-group 
method. This powerful numerical method seems however 
to be limited to inhomogeneous lattice systems at zero 
temperature. In Ref. Owe have pointed out a novel aspect 
of the ID spin-charge separation phenomenon, which ap- 
pears only at finite temperature. Namely, we have shown 
that the propagation of a spin-density packet becomes 
intrinsically diffusive due to the spin drag effect (6] [3 18] 
- the transfer of momentum between fermions of oppo- 
site spin orientations. We have developed a macroscopic 
approach, formulated entirely in terms of collective vari- 
ables, for calculating the propagation of spin-density or 
just density packets, and we believe that the difference 
between the two cases can be experimentally detected. 

However, some difficulties persist, which limit our abil- 
ity to calculate theoretically the dynamics of spin and 
density packets: (i) single-particle thermal excitations, 
which may affect the shape of a spin-density packet in 



ways that are difficult to distinguish from the effect of 
the spin drag, are beyond the reach of the macroscopic 
treatment of Ref. [5] (ii) the macroscopic theory, be- 
ing restricted to long wavelengths and low frequencies, 
misses the interaction of the wavepacket with microscopic 
Friedel-like oscillations, which are induced, for example, 
by the sharp ends of the waveguide; and, finally, (iii) the 
macroscopic theory is restricted to the linear-response 
regime (LRR), whereas any experiment that is likely to 
be performed with ultracold atomic gases would start 
with a strong local disturbance [I]. 

It might appear that these difficulties are formidable 
enough to prevent any further progress, but it is not so. 
It turns out that a theoretical tool already exists, which 
can take care of all three problems in a relatively simple 
manner. It is called "time- dependent spin-current den- 
sity functional theory" (TD-SCDFT) ® EJ1 Q3] and it is 
based on the idea of mapping the time-dependent many- 
body problem into an effective single-particle problem, 
with an effective potential that is designed to produce 
the correct evolution of the spin density and the spin- 
current density. 

Single-particle aspects are naturally included in TD- 
SCDFT because the spin-resolved densities are repre- 
sented as a sum of contributions from one-particle wave- 
functions, which are fully microscopic objects that are 
populated according to a Fermi-Dirac distribution func- 
tion (so thermal effects and Friedel oscillations are auto- 
matically taken into account). Collective effects are in- 
cluded through an exchange-correlation (xc) field, which 
is expressed in terms of collective variables - the spin 
density and the spin-current density - and is derived 
non-empirically from a suitable homogeneous reference 
system. Finally, the theory is not restricted to the LRR. 

In this paper we present the first application of TD- 
SCDFT to the study of the collective dynamics of 
wavepackets in ID Fermi gases, and we compare the 
results with those of Ref. 03 We show that the sim- 
plest xc potential, constructed from the adiabatic local- 



spin-density approximation (ALSDA), is already suffi- 
cient to produce spin-charge separation. But we also ob- 
serve that the density packet gets progressively "contam- 
inated" with microscopic Friedel-likc oscillations coming 
from the ends of the waveguide, which were absent in 
the macroscopic calculation. Going beyond the ALSDA, 
we show that the spin-drag effect is nicely captured and 
continues to be the dominant cause of diffusion. In short, 
we demonstrate the feasibility of a novel and fully micro- 
scopic computational approach, which leads to detailed 
predictions to be compared with future experiments. 
The model — We consider a two-component Fermi gas 
(e.g. a mixture of two hyperfine states of 6 Li [12]) with N 
atoms confined inside a tight waveguide along the x direc- 
tion. The two species of atoms are assumed to have the 
same mass m and different "spin" (hyperfine state label) 
a =| or J,. The fermions interact via a zero-range s-wave 
repulsive potential v(x) = giT>5(x) [13], whose strength 
can be tuned e.g. by using a magnetic field-induced Fes- 
hbach resonance between the two different spin states [2] . 
The Fermi gas is confined in the waveguide by a static 
trap V(x). 

Similarly to what done in Refs. H] and HH we now 
imagine that a local external potential, which couples 
only to atoms of one spin species (hereby referred to 
as "the accumulator"), acts on the Fermi gas for times 
t < creating an accumulation of, say, up-spin atoms 
close to the trap center. The state of the system for 
times t < is denoted by \^o), which is the ground 
state of the Hamiltonian 7i = fi-y + V + V acc - Here 
^y = J2tPl/( 2m ) + ffiD Ei<j <5(*i - Xj) is the Yang 
Hamiltonian [15], V = ^2n V(xi) describes the trap- 
ping potential, and, finally, V acc = J2i,a W <?&) witn 
W a (x) = A a exp[—x 2 /(2w 2 )], is the potential that de- 
scribes the accumulator [16 (Af < and Aj =0). We 
do not make any assumption on the strength A^ of the 
accumulator (i.e. we do not assume to be in the LRR). 

At time t = + the accumulator is suddenly turned off 
and the "charge", n(x,t) — n-\(x,t)+ni(x,t), and "spin", 
s(x, t) — rii(x,t) — n.| (x, t) , densities are allowed to prop- 
agate along the waveguide in the presence of V(x). Here 
n a (x,t) — ( 1 $(t)\ r ipl(x)^p a (x)\ 1 ^(t)) are the spin-resolved 
densities, with |\&(t)) the state of the system at time t [17] 
and ipt-( x ) a field-operator that creates a fermion with 
spin a at position x. 

Dynamics from TD-SCDFT — According to TD- 
SCDFT E3 HI], the spin-resolved densities n a (x,t) 
and the associated current densities j a (x,t) at times 
t > can be found by solving the time-dependent 
Kohn-Sham (KS) equations 
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Vxc [n a , jcr]{x, t) is the KS potential, which includes the 
trapping potential V, the Hartree mean-field potential 
[Vft '[na](x, t) — gi^,n s (x,t), where a = —&], and the xc 
potential - the latter a functional of the spin and cur- 
rent densities. Although, in general, the xc effects in 
TD-SCDFT are represented by a vector potential |10j . 
it turns out that in ID a vector potential can be trans- 
formed into a scalar potential by an appropriate gauge 
transformation: here we have already taken advantage of 
this possibility. The densities are self-consistently deter- 
mined via the usual relation 
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where e a . a are the static KS energies of the initial state 
with chemical potential //. These energies are found by 
solving a static KS self-consistent problem corresponding 
to 7io- The current densities are related to the densi- 
ties by the continuity equations d x j a (x, t) = —dtn a (x, t). 
Note that due to the time-dependence of n a and j a , the 
KS Hamiltonian is time- dependent. 

In order to proceed we need a sensible approximation 
for the xc potential. In order to grasp the physically rel- 
evant facts that occur after the quench described above, 
we make use of the following approximate expression: 
I/ x ( c CT) [n CT , Jff ](x,t) ~ V|lsdaKKM) + V^l[n a , 3a ](x,t). 

In this equation IalsdaKK 3 ''*) represents the ALSDA 
contribution, which depends only on the densities [18.: 
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where £x° m is the xc energy (per particle) correspond- 
ing to the Yang model Tiy, which can be easily found 
from the Bethe- Ansatz equations in the thermodynamic 
limit PUCES]. 

The dynamical (or non-adiabatic) contribution to the 
xc potential, V^l[n a ,j<r]{x,t), is given by 



v£l(x,t) = - 



M) 



dx>F£>(x',t) , 



(4) 



where F^ is the spin-drag-related force [B] exerted by 
the atoms with spin a on the atoms with spin a, 
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n a — *n <r {x,t) 



Due to Galileian invariance this force depends on the rel- 
ative velocity between the two atom species, v a — v s — 
jo/n<7 — js/n s . In Eq. d5J) r s d is the spin-drag relaxation 
time - the inverse of the rate of momentum transfer be- 
tween atoms of opposite spin orientation - which has re- 
cently been calculated in ID [20 . We will use the results 



of Ref. 20 as input for our numerical calculations. Using 
Eq. (I5J) and the continuity equation in Eq. Q we find 
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where T a '{x' \t) — f* dx"d t n a /(x" , t). 

Qualitative analysis of spin-charge separation — Before 
proceeding with the numerical analysis we want to clar- 
ify the mechanism by which the KS equations (TXT) pro- 
duce independent evolutions of the charge and spin den- 
sity. To see this, it is not even necessary to go beyond 
the ALSDA. The essential point is that the KS equa- 
tion guarantees not only the continuity equation but also 
the continuity equation for the momentum density, which 
reads dtjo-(x,t) = —m~ 1 d x P a {x 1 t), where the quantum 
pressure P a (x,t) can be expressed in terms of KS or- 
bitals and is therefore an implicit functional of the den- 
sities. Combining the two conservation laws we arrive 
at dfn a (x,t) — m~ 1 d 2 P (7 (x,t), which looks almost like 
a classical wave equation. Indeed a classical wave equa- 
tion is immediately obtained in the LRR (small deviation 
from homogeneous, unpolarized state) since the quantum 
pressure can then be approximated, in ALSDA, as a lin- 
ear functional of the densities: P a = J2 a > faa'$ n cr' where 
dn a r are deviations from equilibrium. Then, after simple 
algebraic transformation we arrive at two independent 
wave equations for n(x,t) and s(x,t) with two different 
velocities, u n (s) = \A/h ^ hi)/ 171 ' respectively. Admit- 
tedly this analysis pertains to the LLR. Spin and charge 
are not expected to be truly independent in the nonlin- 
ear regime. But since the correct linear-response limit 
is built into the KS equation we expect that a strong 
signature of spin-charge separation will be seen also in 
the nonlinear regime. Our numerical calculations con- 
firm this expectation. 

Numerical results and discussion — We have solved 
Eqs. ([TJ-Q with a two-step predictor-corrector Crank- 
Nicholson scheme. While the scheme outlined above is 
completely general, for simplicity, in the numerical cal- 
culations we have taken a simple box-shaped trapping po- 
tential, V(x) = for -L/2 < x < L/2 and V(x) = +oo 
elsewhere. We use L* = L/N as unit of length and 
E* = h 2 /(mL* 2 ) as unit of energy. For definiteness, we 
fix N T = N i = 40 (N = 80 fcrmions in total), A T = -E*, 
w = L*/V2, g 1D = 2.0 E*L*, and T = 0.5 E*/k B . In 
Fig.[T]wc compare the evolution calculated in the ALSDA 
(which does not include spin drag) with the evolution 
obtained from the macroscopic theory [5]. Both theo- 
ries predict a splitting of the initial peak into two peaks 
that propagate with different velocities (spin being slower 
than charge) in agreement with the general picture of 
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FIG. 1: (Color online) Top panel: charge density n(x, t) (in 
units of 1/L*) as a function of distance x (in units of L*), after 
the subtraction of the background constant rib = N/L. The 
initial density profile is the wavepacket centered at x = 
and labeled by "i = 0". The (red) solid line labeled by 
"MACRO" represents the corresponding n(x, t) at a later time 
t — 5 h/E* , calculated according to the macroscopic theory 
of Ref. H The (blue) solid line labeled by "MICRO" rep- 
resents n(x, t) at the same time but calculated according to 
TD-SCDFT [Eqs. ([T|)-([2j] without spin-drag contribution [i.e. 
with V^l(x,t) = 0]. Bottom panel: spin density s(x,t) (in 
units of 1/L*) as a function of x/L* . Color coding and label- 
ing are the same as in the top panel. 



spin-charge separation. However, in the microscopic cal- 
culation the magnitude of the peaks decreases far more 
rapidly with time. It must be appreciated that this hap- 
pens in spite of the fact that the total spin and parti- 
cle number are conserved quantities in both approaches. 
Note also that the width of the microscopic results in 
Fig. [T] is much larger than that of the macroscopic re- 
sults: the reason is that the microscopic theory includes 
diffusive-like mechanisms related to single-particle ther- 
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FIG. 2: (Color online) Same as in the bottom panel of Fig. n] 
but with spin-drag contribution. In the inset we show a 3D 
plot of the TD-SCDFT result for the spin density s(x,t) (in 
units of 1/L*) as a function of x/L* and tE* /h. 




FIG. 3: (Color online) A 3D plot of the TD-SCDFT result 
for the charge density n(x,t) (in units of 1/L*) as a function 
of x/L* and tE* /h. 



mal excitations [see Eq. (pi], which operate both in the 
charge and spin channels and are completely missed by 
the macroscopic theory. Finally, notice the asymmet- 
ric forward- leaning shape of the density pulse calculated 
from the microscopic theory. This is a nonlinear effect, 
likely due to the fact that the local velocity, proportional 
to the density, is higher at the center of the pulse than 
at its edges [2"T] . 



Fig. [2] shows the effect of spin drag, which enters the 
macroscopic calculation as a dissipative term in the wave 
equation for s(x,t) [see Eq. (10) in Ref. |5], the micro- 
scopic one through the current-dependent part of the xc 
potential [Eq. |6])]. It is evident that the effect of the spin 
drag is much more pronounced in the microscopic cal- 
culation, where the double-peak structure is completely 
lost. This happens in spite of the fact that the spin-drag 
coefficient r s( j is generally smaller in the microscopic cal- 
culation, due to the effect of the spin polarization |20j . 
Clearly, the inclusion of microscopic excitations weakens 
the collective behavior of the spin density. The large 
difference between the macroscopic and microscopic re- 
sults demonstrates the importance of relying on the latter 
when performing quantitative calculations. 

In Fig. [3] we show a 3D plot of the time evolution of 
the density packet. What is notable here is that already 
at short times some density waves appear, coming from 
the sharp edges of the waveguide: at larger time they 
mix with the original packet. These Friedel-like oscilla- 
tions are completely beyond the power of the macroscopic 
approach. In conclusion, the above calculations amply 
demonstrate the versatility of the TD-SCDFT method in 
producing detailed results which can be compared with 
future experiments on the propagation of charge and spin 
pulses in ultracold Fermi gases. 
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